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Abstract. We provide microscopic diagrammatic derivations of the Molecular Co- 
herent Potential Approximation (MCA) and Dynamical Cluster Approximation 
(DCA) and show that both are ^-derivable. The MCA (DCA) maps the lattice 
onto a self-consistently embedded cluster with open (periodic) boundary condi- 
tions, and therefore violates (preserves) the translational symmetry of the original 
lattice. As a consequence of the boundary conditions, the MCA (DCA) converges 
slowly (quickly) with corrections 0(1/L C ) (C(l/L 2 )), where L c is the linear size 
of the cluster. However, local quantities, when measured in the center of the MCA 
cluster, converge more quickly than the DCA result. These results are demonstrated 
numerically for the one-dimensional symmetric Falicov-Kimball model. 



Introduction 

One of the most active areas in condensed matter physics is the search for 
new methods to treat chemically disordered and correlated systems. In these 
systems, especially in three dimensions or higher, approximations which ne- 
glect long ranged correlations are generally thought to provide a reasonable 
first approximation for many properties. 

Perhaps the most successful of these methods are the Coherent Poten- 
tial Approximation (CPA) jD and the Dynamical Mean Field Approxima- 
tion (DMFA) for disordered and correlated systems, respectively. 
Although these approximations have different origins, they are formally re- 
lated. Both are single site theories where non-local effects of disorder and 
correlations are treated in a mean field approximation. Diagrammatically, 
both the DMFA |] and the CPA § may be defined as theories which com- 
pletely neglect momentum conservation at all internal vertices. When this 
principle is applied, the diagrammatic expansion for the irreducible quanti- 
ties in each approximation collapses onto that of a self-consistently embedded 
impurity problem. 

Many researchers have actively searched for a technique to restore non- 
local corrections to these approaches. Here, we discuss just two approaches 
which are fully causal and self-consistent: the Molecular Coherent Poten- 
tial Approximation (MCA) [^J|] and the Dynamical Cluster Approximation 
(DCA) P,p]P,|§. Recently the Cellular Dynamical Mean Field Approach 
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[ fL2| was proposed for ordered correlated systems, while the Molecular Co- 
herent Potential Approximation has traditionally been applied to disordered 
systems. Since both methods share a common microscopic definition we use 
the term MCA to refer to both techniques in the following. 

While the MCA is traditionally defined in the real space of the lattice, 
the DCA is traditionally defined in its reciprocal space. In the MCA, the sys- 
tem lattice is split into a series of identical molecules. Interactions between 
the molecules are treated in a mean-field approximation, while interactions 
within the molecule are explicitly accounted for. In the DCA, the reciprocal 
space of the lattice is split into cells, and momentum conservation is neglected 
for momentum transfers within each cell while it is (partially) conserved for 
transfers between the cells. These approximations share many features in 
common: they both map the lattice problem onto that of a self-consistently 
embedded cluster problem. Both recover the single site approximation (CPA 
or DMFA) when the cluster size reduces to one and become exact as the clus- 
ter size diverges. Both are fully causal and provided that the clusters 
are chosen correctly || , they maintain the point group symmetry of the orig- 
inal lattice problem. Here, we provide a microscopic diagrammatic derivation 
of both the MCA and the DCA, and explore their convergence with increasing 
cluster size. 

Formalism 

We will employ a diagrammatic formalism to derive the MCA and DCA, 
assuming that a collection of electrons on a lattice, with Green function G(k) 
interact through an interaction V(k). 

From the Lattice to the Cluster 

Since our object is to define cluster methods, we divide the original lattice 
of N sites into N/N c clusters (molecules), each composed of iV c = sites, 
where D is the dimensionality. L c need not be an integer, for example in 
a two-dimensional square lattice, a diamond cluster with N c = 8 will have 
L c = 2\2. However, care must be taken so that the clusters preserve the point 
group symmetry of the original lattice. We use the coordinate x to label the 
origin of the clusters and X to label the N c sites within a cluster, so that 
the site indices of the original lattice x = X + x. The points x form a lattice 
with a reciprocal space labeled by k. The reciprocal space corresponding to 
the sites X within a cluster shall be labeled K, with K a — n a ■ 2ir/L c and 
integer n a . Then k = K + k. Note that e = 1 since a component of x 
must take the form m a L c with integer m a . 

The formal mapping between the full lattice problem and the cluster is 
accomplished by relaxing the condition of momentum conservation at the 
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Fig. 1. Definitions of the parameters in real (left) and reciprocal (right) space. 
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internal vertices of the compact diagrams. Momentum conservation at each 
vertex is described by the Laue function 

A = e ix.(k 1+ k 3+ ...,-ki-i4-...) = ^ kl+k2+ ..., k , +k , + ... , (l) 

X 

where ki, k2 (k' 1; k' 2 ) arc the momenta entering (leaving) the vertex. Miillcr- 
Hartmann (3) showed that the Dynamical Mean Field (DMF) theory may 
be derived by completely ignoring momentum conservation at each internal 
vertex by setting A = 1. Then, one may freely sum over all of the internal 
momentum labels, and the graphs for the generating functional <f> and its 
irreducible derivatives, contain only local propagators and interactions. 

The DCA and MCA techniques may also be defined by their respective 
Laue functions. In the MCA, we approximate the Laue function by 

A M c = e 4X - (Kl+ltl+K2+t2+ -~ K 'i^'i" K 2- lt 2--) . (2) 
x 

Thus, the MCA omits the phase factors e lk x resulting from the position of 
the cluster in the original lattice but retains the (far less important) phase 
factors e lk x associated with the position within a cluster. In the DCA we 
also omit the phase factors e lk x , so that 

Adc = N c 5-k. 1 +k 2 +- i k.{+tc' 2 +- • (3) 

Both the MCA and DCA Laue functions recover the exact result when N c — > 
00 and the DMFA result, A = 1, when N c = 1. 

If we apply the MCA Laue function Eq. ^| to diagrams in <P, assuming 
that V is a two-particle interaction then each Green function leg is replaced 
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by the MCA coarse-grained Green function (we have dropped the frequency 
dependence for notational convenience) 

G(X 1 ,X 2 ;x = 0) = 

Kl,K 2 

ki,k 2 

N 2 

j± GpLuXiMM), (4) 

ki,k 2 

or in matrix notation for the cluster sites Xi and X2 

G=^£6(k), (5) 

k 

since G can be chosen diagonal in ki , k2 . Similarly each interaction line is 
replaced by 

K(X 1 ,X 2 ;x = 0) = 

_L_£ e^+kO-x^^^^.^^^^CK.+k^.x. = 

Kl,K 2 

ki,k 2 

N 2 

j± nx^Xa.kx.ka), (6) 

ki,k 2 

or in matrix notation for the cluster sites Xi and X2 

v = (7) 
ic 

The summations of the cluster sites X within each diagram remain to be 
performed. G and V are propagators which are truncated outside the cluster. 
I.e., if the interaction V is non-local, V will include only interactions within, 
but not between, clusters. Thus the inclusion of the phase factors e lk ' x in the 
MCA Laue-function Eq. || leads directly to a cluster approach formulated in 
real space that violates translational invariance. Therefore the Green function 
and interaction are functions of two cluster momenta Ki, K2 or two sites Xi, 
X2 respectively. 

If we apply the DCA Laue function Eq. g Green function legs in <P are 
replaced by the DCA coarse grained Green function 

k 
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since Green functions can be freely summed over the k vectors within a cell 
about the cluster momentum K. Similarly, the interactions are replaced by 
the DCA coarse grained interaction 

V(K) = |^F(K,k). (9) 

k 

As with the MCA, the effect of coarse-graining the interaction is to reduce 
the effect of non-local interactions to within the cluster. The resulting com- 
pact graphs are functionals of the coarse grained Green function G(K) and 
interaction V(K), and thus depend on the cluster momenta K only. For ex- 
ample, when N c = 1 , only the local part of the interaction survives the coarse 
graining. As with the MCA, within the DCA it is important that both the 
interaction and the Green function are coarse-grained |T(|. In calculations 
where a non-local interaction is not coarse-grained, poor results are obtained 




Prom the Cluster to the Lattice 

To establish a connection between the cluster and the lattice we minimize 
the lattice free energy 

F = -fc B T(<5-tr[EG] +trln[G]) (10) 

where is the generating functional composed of all closed compact (single- 
particle irreducible) graphs, £ is the lattice self-energy and G is the full 
lattice Green function. The trace indicates summation over frequency, mo- 
mentum and spin. As discussed in many-body texts 0], the additional free 
energy due to an interaction may be described by a sum over all closed con- 
nected graphs. These graphs may be further separated into compact and 
non-compact graphs. The compact graphs, which comprise the generating 
functional consist of the sum over all skeletal graphs (those with no in- 
ternal parts representing corrections to the single-particle Green function). 
The remaining graphs comprise the non-compact part of the free energy. In 
the infinite-dimensional limit, <P consists of only local graphs, with non-local 
corrections of order l/D. However, for the non-compact parts of the free en- 
ergy, non-local corrections over arbitrary lengths are of order one, so the local 
approximation applies only to 

To see this, consider the simplest non-local corrections to non-compact 
and compact parts of the free energy of a Hubbard-like model, illustrated in 
Fig. |^. Here the upper (lower) circle is a set of graphs composed of intrasite 
propagators restricted to site n (the origin). Consider all such non-local cor- 
rections on the shell of sites which are n mutually orthogonal unit translations 
from the origin. In the limit of high dimensions, there are 2 n D\/((D— n)!n!) ~ 
0(D n ) such sites. Since as D — > oo, G(r) ~ D~ T I 2 [0, the legs on the compact 
correction contribute a factor 0(D~ 2n ) whereas those on the non-compact 
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Fig. 2. Non-compact (left) and compact (right) non-local corrections to the free 
energy functional. Here the upper (lower) circle is meant to represent a set of 
graphs which are closed except for the external lines shown, and restricted to site 
n (the origin). 



correction contribute 0{D~ n ), Therefore the compact non-local correction 
falls as D~ n and is very short-ranged; whereas, the non-compact correction 
remains of order one, regardless of how far site n is from the origin Jig ]. As 
we will see below, the essential approximation of the DCA and the MCA is 
to use the cluster propagators, which are accurate only for short distances, 
to construct various diagrammatic insertions. In high dimensions, or in finite 
dimensions when the Green functions fall exponentially with distance, this is 
a good approximation for the compact graphs which comprise ^, but a poor 
approximation for the non-compact graphs []l5f . 

Thus, we will approximate the generating functional <P with its cluster 
counterpart <P C by replacing the Laue function with either A uc or A mc , but 
this approximation will not be used in the parts of the free energy coming 
from non-compact graphs. The free energy then reads 

F = -k B T{<P c - tr[£G] +trln[G]) (11) 

F is stationary with respect to G when |^ = 0. This happens for the MCA 
if we estimate the lattice self energy as 

J- e - 4 ( Kl +^)- x ^ MC .(X 1 ,X 2 )e 4 ( K2+ ^)- x 2. (12) 

Xi ,X2 

Thus, the corresponding lattice single-particle propagator reads in matrix 
notation 

G(k, z) = [zl - e(k) - Smc(z)] 1 , (13) 

where the dispersion e(k) and self-energy Emc(z) are matrices in cluster real 
space with 

[e(k)] Xl x 2 =e(X 1 -X 2 ,k) (14) 
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being the intracluster Fourier transform of the dispersion. For the DCA, 
-S'(k) = Z , /3C'(K) is the proper approximation for the lattice self energy cor- 
responding to <1>dc- The corresponding lattice single-particle propagator is 
then given by 



G(K,k;z) 



1 



- e K+k - Zdc{K,z) 



(15) 



Both the MCA and DCA are optimized when we equate the lattice and cluster 
self energies. A similar relation holds for two-particle quantities. Thus, with 
few exceptions , only the irreducible quantities on the cluster and lattice 
correspond one-to-one. 

The MCA (DCA) algorithm, illustrated in Fig. |, follows directly: We 
first make an initial guess for the cluster self-energy matrix S. This is used 
with Eqs. || and [l3| (|| and |l5|) to calculate the coarse-grained Green function 

G. The cluster excluded Green function G = [G + Smc\~ X (Go(K) = 
[G(K) _1 + Z£>c(K)] _1 ) is defined to avoid overcounting self energy correc- 
tions on the cluster. It is used to compute a new estimate for the cluster 
self-energy which is used to reinitialize the process. 



MCA Algorithm 



_ Cluster 
Solver 




GJK 




GJK) 



--Z(K) 



DCA Algorithm 

Cluster 
Solver 



G(K) 




G(G(k)) 

Fig. 3. The MCA and DCA algorithms. The two differ mainly in that in the DCA, 
the Green functions are diagonal in K, while in the MCA, they are matrices in the 
cluster coordinates X. 



Once convergence is reached, the irreducible quantities on the cluster may 
be used to calculate the corresponding lattice quantities. For example, the 
cluster self energy and irreducible vertex functions may be used in the Dyson 
and Bethe-Salpeter equations to calculate the Green function and suscep- 
tibilities. In order to obtain a smooth DCA self energy for the calculation 
of the band structure and spectra, it may be interpolated into the Brillouin 
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zone of the lattice (however, such interpolation should be avoided during the 
self-consistency loop, as it may lead to causality violations). The optimal 
MCA self energy, Eq. [l2] is a function of two momenta since the translational 
invariance of the lattice is violated. Kotliar et al., have introduced a cluster 
averaging scheme to obtain a self energy as a function of one momenta which 
may be employed after convergence is obtained [12] . 



The Small Parameter, P: the Coupling Between the Cluster and 
its Host 

In order to compare the character of the two different cluster approaches as 
a function of the cluster size N c it is instructive to rewrite the corresponding 
coarse grained Green-functions Eqs. || and || to suitable forms by making use 
of the independence of the self-energy S on the integration variable k. For 
the MCA coarse grained Green function we find 

- r ~ i _1 

G(z) = zl-e a - Smc(z) - r MC (z) , (16) 

with the "cluster-local" energy i a = N c /N^2^e(li). ^ or ^ ne DCA we obtain 
a similar expression 

G(K, z) = [z- e K - S DC {K, z) - r DC (K, z)]" 1 , (17) 

with the coarse grained average ck = N c /N e (K, k). The hybridization 
functions Pmc / dc{ z ) describe the coupling of the cluster to the mean- field 
representing the remainder of the system. 

The behavior of r for large N c is important. For the MCA, r averaged 
over the cluster sites and frequency 



Xi ,x 2 



(18) 



where L c = N^ D is the linear cluster size. A detailed derivation of this 
form is presented in the appendix. However, since in the MCA the cluster 
is defined in real space with open boundary conditions, this form is evident 
since only the sites on the surface cx 2D ■ Lf^ 1 of the cluster couple to 
the effective medium and N c = . For the DC A we show in the appendix 
that r(K) ~ 0(l/Nc^ D ) (see also |ll[) so that we obtain for the average 
hybridization of the DCA cluster to the effective medium 

r °c = W c ^ FDc{K) ~°{ll) ■ (19) 

The DCA coarse graining results in a cluster in K-space; thus, the corre- 
sponding real space cluster has periodic boundary conditions, and each site 
in the cluster has the same hybridization strength _T with the host. 
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As shown in the appendix, the average hybridization strength r acts as 
the small parameter in both the MCA and the DCA. Thus the MCA (DCA) 
is an approximation with corrections of order P ~ 0(1/L C ) (~ 0(1/L%)). 



Numerical Results 

To illustrate the differences in convergence with cluster size N c we performed 
MCA and DCA simulations for the symmetric one-dimensional (ID) Falicov- 
Kimball model (FKM). At half filling the FKM Hamiltonian reads 

H = -t EKH+i + h - c -) + Uj2( n >- V2)(nf - 1/2) , (20) 

i i 

with the number operators nf = d\d i and n{ = fjf i and the Coulomb 
repulsion U between d and / electrons residing on the same site. The FKM 
can be considered as a simplified Hubbard model with only one spin-species 
(d) being allowed to hop. However it still shows a complex phase diagram 
including a Mott gap for large U and half filling, an Ising-like charge ordering 
with the corresponding transition temperature T c being zero in ID, and phase 
separation in all dimensions. The bare dispersion (in ID) ej. = 2tcosk; thus 
for t = 1/4 the bandwidth W — 1 which we use as unit of energy. To simulate 
the effective cluster models of the MCA an the DCA we use a quantum Monte 
Carlo (QMC) approach described in Jl0| . 

To check the scaling relations Eqs. [18[and [l9| we show in FigJJ the average 
hybridization functions Pmc and Pdc for the MCA and DCA respectively at 
the inverse temperature (3 = 17 for U = W = 1. For N c = 1 both approaches 
are equivalent to the DMFA and thus Pmc — Pdc- For increasing N c Pmc 
can be fitted by 0.3361/iV c and P DC by 1.1946/ATf when N c > 2. Cluster 
quantities, such as the self energy and cluster susceptibilities, are expected 
to converge with increasing iV c like P. This is illustrated in the inset for the 
staggered (Q = n) charge susceptibility Xc(Q) of the cluster. 

Since only the compact parts represented by <P of the lattice free energy 



(Eq. fO) are coarse-grained, this scaling is expected to break down when lat- 
tice quantities, such as the lattice charge susceptibility, are calculated. The 
susceptibility of the cluster Xc{Q) cannot diverge for any finite A^ c ; whereas 
the lattice x(Q) diverges at the transition temperature T c to the charge or- 
dered phase. Note that the residual mean-field character of both methods can 
result in finite transition temperatures T c > for finite N c < oo. However as 
N c increases, this residual mean field character decreases gradually and thus 
increased fluctuations should drive the solution to the exact result T c = 0. 

In the DCA M], xiQ) is calculated by first extracting the corresponding 
vertex function from the cluster simulation. This is then used in a Bethe- 
Salpeter equation to calculate x(Q)- T c is calculated by extrapolating x(<3) -1 
to zero using the function x(Q) _1 oc (T— T c ) 7 (see inset to Fig.^J) . This proce- 
dure is difficult, if not impossible, in the MCA due to the lack of translational 
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Fig. 4. The average integrated hybridization strengths t of the MCA (squares) 
and DCA (circles) versus the cluster size N c when (3 = 17 and U = W = 1 in 
the symmetric model. The solid and dashed lines represent the fits 1.1946/JV^ and 
0.3361/iV c respectively. Inset: Convergence of the cluster charge susceptibility for 
Q — tv. The solid and dashed lines are quadratic and linear fits, respectively. 

invariance. Here, we calculate the order parameter m(T) — 1 /N c 1)* ( n i ) 
in the symmetry broken phase. T c is then obtained from extrapolating m(T) 
to zero using the function m(T) cx (T c — T) 13 . For the DCA this extrapolation 
is shown by the solid line in the inset to Fig.|] for N c — 4. The values for 
T c obtained from the calculation in the symmetry broken phase and in the 
unbroken phase must agree, since as we have shown above, both the DCA 
and MCA are ^-derivable. This is illustrated in Fig.|| for the DCA. 

A comparison of the DCA and MCA estimate of T c is presented in Fig. |[ 
T c obtained from MCA (squares) is larger than T c obtained from DCA (cir- 
cles). Moreover we find that the DCA result seems to scale to zero almost 
linearly in 1/N C (for large enough N c ), whereas the MCA does not show any 
scaling form and in fact seems to tend to a finite value for T c as N c — > oo. 
This striking difference of the two methods can be attributed to the differ- 
ent boundary conditions. The open boundary conditions of the MCA cluster 
result in a large surface contribution so that -Tmc > ^~!dc- This engenders 
pronounced mean field behavior that stabilizes the finite temperature transi- 
tion for the cluster sizes treated here. For larger clusters we expect the bulk 
contribution to the MCA free energy to dominate so that T c should fall to 
zero. 

Complementary results are found in simulations of finite-sized systems. 
In general, systems with open boundary conditions are expected to have a 
surface contribution in the free energy of order 0(1/L C ) fl7| . This term is 
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Fig. 5. The transition temperature T c for the DCA (circles) and MCA (squares) 
when U = W = 1 versus the cluster size N c . For all values of N c the DCA prediction 
is closer to the exact result (T c = 0). Inset: Order parameter m(T) and inverse 
charge susceptibility x(Q) _1 versus temperature. The solid (dashed) line represents 
a fit to the functions m(T) ex (T c - T) with j3 = 0.245 (x{T) oc (T — T c )~ 7 with 
7 = 1.07). 



absent in systems with periodic boundary conditions. As a result, simulations 
of finite-sized systems with periodic boundary conditions converge much more 



quickly than those with open boundary conditions 18 . 

Thus far, we have shown that the DCA converges more quickly than the 
MCA for critical properties and for extended cluster quantities (e.g. the clus- 
ter susceptibility). This is due to the differences in the boundary conditions, 
and the coupling to the mean-field. Whereas each site in the DCA experi- 
ences the same coupling to the mean- field host, in the MCA only the sites 
on the boundary of the cluster couple to the host. Provided that the system 
is far from a transition, the sites in the center of the cluster couple to the 
mean-field only through propagators which fall exponentially with distance. 
Thus, one might expect that local results, such as the single-particle density 
of states, might converge more quickly within the MCA provided that they 
are measured on these central sites. This is illustrated in Fig. ^ where we plot 
the single-particle density of states calculated with the DCA and the MCA 
on the two central sites. 



Summary. 



By defining appropriate Laue functions, we provide microscopic diagram- 
matic derivations of the MCA and DCA. We show that they are <£-derivable, 
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Fig. 6. Comparison of the single-particle density of states (DOS) calculated with 
the DCA (left) and the MCA (right) on the two central sites for various cluster 
sizes N c when j3 = 10 and U = W = 1 (since the DOS is symmetric around u = 0, 
only the to > part is shown). Local quantities, such as the DOS, converge more 
quickly when calculated on the central sites of the MCA cluster than the cluster 
averaged DCA result. 

and that the lattice free energy is optimized by equating the irreducible quan- 
tities on the lattice to those on the cluster. The MCA maps the lattice to a 
cluster with open boundaries and consequently, the cluster violates transla- 
tional invariance. In contrast, the DCA cluster has periodic boundary condi- 
tions, and therefore preserves the translational invariance of the lattice. 

These differences in the boundary conditions translate directly to different 
asymptotic behaviors for large clusters N c . As we find analytically as well 
as numerically, the surface contributions in the MCA lead to an average 
hybridization F of the cluster to the mean field that scales like l/L c as 
compared to the scaling of the DCA. Since F acts as the small parameter 
for these approximation schemes, the DCA converges much more quickly than 
the MCA. These effects are more pronounced near a transition, where the 
large surface contribution of the MCA stabilizes the mean-field character of 
the transition. Consequently, the DCA result for the transition temperature 
T c of the ID symmetric FKM model scales almost like 1/N C to the exact 
result T c = 0, whereas the MCA result converges very slowly. The boundary 
conditions also differ in that only the MCA sites at the surface of the cluster 
couple to the mean field; whereas, all DCA cluster sites have an equal coupling 
to the mean field host. As a result local quantities, such as the density of 
states, when measured on the central sites of the MCA cluster converge more 
quickly than than corresponding DCA results. 

Thus, for critical properties and extended cluster quantities, the DCA 
converges far more quickly than the MCA; whereas for local quantities which 
may be measured at the central sites of the MCA cluster, the MCA converges 
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more quickly. Since the origin of these differences lies in the different boundary 
conditions we expect them to hold generally for any model of electrons moving 
on a lattice. 
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Appendix 

To differentiate the MCA from the DCA we find expressions for the cor- 
responding host-functions r. To this end we split the hopping integral t 
into an inter- and intracluster part. For the MCA the intercluster hopping 
Pmc]xi,x 2 = 7mc(Xi,X 2 ) is defined as T M c{^) = e(k) - e Q and reads in 
real space 

f MC (x) = f(x)-e 5(x). (21) 

Since e G = f(x = 0), the intercluster hopping matrix Tj;c(x) is only finite 
for x ^ as expected. Thus, Tmc(x) has non- vanishing matrix-elements 
only for sites on the boundary of the cluster. For the DCA the intercluster 
hopping integral is analogously defined as Tbc(K, k) = e(K, k) — In and can 
be written in real space as 

Tdc(x) = t(x) - t«S c (x) , (22) 

where [i]xj.x a = *(Xi-X 2 ) = 1/N C J2 K e iK(Xl_X2) e K and the cluster delta- 
like function [<5 c (x)] Xl x 2 = <5 c (Xi - X 2 ,x) = N c /N e *-(Xi-X 2 +x) , We 
have discussed elsewhere that 5 c (Xi — X 2 ,x) ~ 1 for x = and falls off 
rapidly for finite x. It is important to note that t(Xi — X 2 ) = t(Xj — X 2 , x = 
0) for N c = N, but |t(Xi - X 2 )| < |i(Xi - X 2 ,x = 0)| for A^ c < N. It follows 
that the DCA intercluster hopping Tdc*(Xi — X 2 ,x) for x = 0, i.e. between 
sites belonging to the same cluster, is even finite when A^ c < A^; however it is 
strongly reduced compared to x ^ 0. This is a consequence of the periodicity 
of the DCA cluster. Therefore we find that the effective model of the DCA 
cannot be thought of N/N c clusters composing the original lattice, but we 
have to rather think of it as a renormalized model in real space. 

Now we proceed with finding expressions for the hybridization functions 
r defined in Eq. @ and @ for the MCA and DCA. We start with the 
MCA. To this end we define a Green function matrix 

g(z) = \zI-e-E MC (zj\ ' , (23) 
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that is localized on the impurity cluster. With this definition Eq. (0) reads 
N 2 

1 y C 

N 2 



« 7U2 r , _1 

^ = ^IE 6 W- g-'-Tuc , (24) 



k 



and Eq. ( |l3| ) can be rewritten to 

G(k) = § + §f M c(k)G(k). (25) 

After inserting Eq. ( p5|) into Eq. (^J) and using X)k^Mc(k) = we obtain 
after some algebraic transformations for the MCA-hybridization matrix 



r, 



MC 



7 +|E f ^(k)G(k) 



N 

k 



-1 

X 



^^f M c(k)G(k)f MC (k). (26) 

k 

Since, as we pointed out above, T has non-zero matrix-elements only between 
sites on the boundary of the impurity cluster, _T is finite only for sites on 
the boundary but vanishes for sites inside the cluster. Thus, the average 
hybridization strength of the MCA cluster per site 

^ = wH A/c(Xx,X 2 ) ~ O ( 2 -f) , (27) 

iVc X!,X 2 ^ LcJ 

since only the sites on the surface oc 2D ■ L®^ 1 of the cluster couple to the 
effective medium and N c = . 

For the DCA we can follow the steps presented above for the MCA since 
-Tdc^K) can be considered a diagonal matrix in K. We obtain for the DCA 
hybridization function 

^5>£ c (K,k)G(K,k) 

r DC (K) = ™* . (28) 

l + ^E T ^ K ' k ) G ( K ' k ) 

k 

It should be stressed that due to the periodicity of the DCA cluster ev- 
ery site of the impurity cluster couples to the effective medium. The mo- 
menta k are restricted to a DCA coarse graining cell and therefore maximal 
of the order 0{Ak), where Ak oc 1/L. Since e(K) - e K ~ 0({Ak) 2 ), we 
find by performing a Taylor series expansion of Tdc^K, k) around e(K) that 
TD C (K,k) ~ 0(1/L C ). For the overall hybridization of the DCA cluster to 
the effective medium we thus obtain 

^ = ^E r ^( K )~°(zf) • ( 29 ) 
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In both the DCA and the MCA, the average hybridization strength acts 
as the small parameter. The approximation performed by the DCA (MCA) is 
to replace the lattice Green function G(K, k, z) = [z — e K , j — S(K, k, z)]^ 1 

(G(k) = [zl — e(k) — £(k, by its coarse grained quantity G(K) (G) in 

diagrams for the generating functional <P. 

According to Eq. ( ^7|) the coarse grained Green function of the DCA can 
be expressed as G(K, 

*) = [z-e K - ^dc(K, z) - r DC {K, z)] \ For the 
time being, assume that £dc(K, z) has corrections of the same, or higher, 
order in 1/ L c as the average hybridization strength. Both the self energy and 
e K+ j c = 6k + 0(1/ L c ) with the leading order corrections being linear in k. 
Since furthermore r DC (K) ~ 0(1/1% G(K, k) a G(K) +0(1/L C ). The 
diagrams for <P however are summed over k, so that the terms ~ 0(1/ L c ) 
coming from e K+ ^ and similar terms from the self energy vanish and only 
the terms ~ 0(1/L^), or higher, survive. Thus we find for the DCA that 
<P « <Pdc + 0(1/ L^). The corresponding estimate of Z , £>c , (K, z) will also 
have corrections 0(1/L?), confirming the assumption above. 

According to Eq. ( fig ) , the coarse grained propagator of the MCA can be 

written as G(z) = [zl — e a — £mc(z) ~ Afc(z)]~ • If we assume that Emc 
converges with L c as fast or faster than 1 / L c we see with the L c -dependence 
of f MC in Eq. @ that G(k) » G + G(l/L c ) since e(k) = e + 0(1/L C ). 
Thus we obtain for the MCA approximation <P rj <S>mc +0(1/L c ) and £ will 
converge with L c as 0(1/ L c ) confirming our assumption. 
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